The Decoherence of the Electron Spin Polarization and 
Meta-stability of ^^C Nuclear Spins in Diamond 

P.R. Crompton 
Department of Applied Maths, School of Mathematics, 
University of Leeds, Leeds, LS2 9JT, UK. 
(Dated: June 26, 2009) 

Abstract 

Following the recent successful manipulation of entangled ^^C atoms on the surface of Diamond 
we calculate the decoherence of the electron spin polarization in Diamond via a nonperturbative 
treatment of the time-dependent Greens function of the Central-Spin model, describing the Hyper- 
fine couplings of the electron to a bath of interacting ^^C atoms, for arbitrary initial polarizations, 
applied field strengths, and for up to eight entangled ^^C atoms. We compare our numerical re- 
sults, which are regularized by ^-function prescription, with the exact treatment available that is 
available in the fully initially polarized limit of the non-Markovian dynamics regime, and comment 
on the relevance of dangerously irrelevant ^^C dipole-dipole interactions to the meta-stability of 
^■^C nuclei flip-flop processes. 
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Recently, several experimental groups [ll-Q] have investigated the feasibility of processing 
Quantum Information via the manipulation of optically excited electron spins |5| in Dia- 
mond. These studies focus on the Diamonds' optically active Nitrogen Vacancy centers and 
the triplet that is formed there through the Zeeman splitting of a Nitrogen impurity coupled 
to a lattice vacancy There are two known charged states of this point defect, which are 
neutral and negatively charged, respectively, with spin states of rris = and = 1. Hence, 
stable qubits have beed formed in Diamond consisting of a single electron spin in the = 1 
state coupled to a bath of the neighbouring nuclear spins which slowly decoheres via the 
spectral diffusion of the nuclear spin polarizations. Diamond offers several advantages over 



conventional Quantum Dot materials, such as GaAs for this purpose j?! 
the density of states of this system is relatively low (with a low electron-phonon coupling), 
whilst the oscillator strength for the electron dipole transition is relatively large. This is true 
in particular for the electron's coupling to the ^^C atoms which have a natural abundance of 
1.1% 0]. Consequently, it has been reported that electron spin-dephasing times of 0.5 — 
can be achieved in Diamond, at 300K - without the need for cryogenic cooling jH-Q], which 
makes Diamond an attractive candidate for Quantum Information Processing. 

This spin decoherence of single electrons in the = 1 state of the Nitrogen Vacancy cen- 
ters of Diamond can thus be modelled via the competition between; the spin-orbit coupling 
of the electron and ^^C nuclei, the Hyperfine coupling between the electron and a spin-bath 
of ^^C nuclear spins, and the dipole-dipole interaction between the ^^C bath nuclei [9l, which 

nn □ 

together is given via the following Central-Spin model Hamiltonian [6||10|-[12l|. 

N N N N 

i7 = e55^ + 5^6,/f + + J2 'f [{Un,,){I,.n,,)-UI,] (1) 

1=1 i=l 1=1 i^j.j=\ V 

where S^{If) is an operator for the longitudinal spin of the electron in the adjoint 
representation (zth nuclear spin), /j is an operator for the ith ^^C nucleus in the spin- 
bath, is the number of bath nuclei and unit vector linking the centers of the 
ith and jth ^^C nuclear dipole moments with separation Vij. Since the relative strength 
of the couplings of the Zeeman, Hyperfine and dipole-dipole interactions in Diamond are 



. This is because 



{es = 115.QfieVT-\ej = 0.1196 fieVT-\ Ai = 0.09873 /leV, e^d = 1.2566 x 10~^ fieV }\li, 
naively, the dipole-dipole ^^C interactions can be neglected. However, what makes this 
model very interesting is the difference between the = and = 1 states in the applied 



2 



field, since for rus = tlie nuclear spins all precess with the same frequency, whilst for 
ms = 1 the bath nuclei spins diffuse via the Hyperfine interaction Hence, although we 
choose to represent the triplet state of the electron via the adjoint representation, if the 
field strength is tuned to the level at which the energies of the = and = 1 states 
become equal (~ 500(7), the otherwise irrelevant ^^C dipole-dipole can become important 
to the decoherence of the electron spin through quantum tunnelling between the triplet and 
singlet state. In such a case the Hyperfine intercation of the electron-dipole interaction 
should then lead to a metastability in the ^^C nuclei spin polarizations. 

Via Sturm-Liouville theory (in the absence of the dipole-dipole interaction), the longi- 
tudinal and transverse spin components of a single electron in the Central-Spin model are 
exactly separable, which yields the following Nakajima-Zwanzig generalized master equa- 
tions (GMEs) 3, 

{S^)t = N,{t)-z f dt'T.,,{t-t'){S'),, (2) 
Jo 

{S+)t = tuj^{S+)t-t f dt'E++{t-t'){S+)t' (3) 

^0 

where Nz, T,zz and S4.+ are given by matrix elements of the reduced self energy S(t — t'), 
ujn = B{es—ej)+hj, with hj an eigenvalue of the initial angular momentum eigenstates of the 
bath, and B is the applied field strength. From which, the time dependence of the electron 
spin polarization can be evaluated via Laplace transform. Importantly, even for this simple 
case, the above GMEs are nonperturbative for Un/N > 1 and so the self-energy contributions 
do not resum and grow as a function of time. To make this GME analytically tractable in 
the pertubative regime (for arbitrary initial polarizations in the high field limit), the above 
self energies can be reparameterised as continuous functions of the Hyperfine couplings (via 
Ak x). In the Born approximation (with S^^ = Yl^^z {I±{s))) this yields the functions, 

I±is) = TT7 > 7T- I±is) = — / dx- —, Z/ = 1, (4) 

k 2 

where Pi{m) is the probability of finding a nuclear spin / with z-projection m, and d 



the dimension of the system [13|]. Hence, for u < 1 the electron spin decoherence envelope 
is oscillatory and is governed by nuclei close to the central electron (Markovian), whereas 
for z/ > 1, the electron spin decoherence envelope is slowly decaying and governed by the 
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nuclei farthest away (non-Markovian). However, in the field splitting problem we consider in 
(1) we need to take a more general approach to evaluating these integrals, since delineating 
the regimes via v is adversely affected by the quantum tunnelling, as the Markovian and 
non-Markovian regimes can become interchanged [l^. In a renormalization group sense, 
z/ is a dangerously irrelevant scaling parameter through the flip-flop dynamics of the ^^C 



nuclei 



15|(l6|[l7|. 



To tackle the fleld splitting case in (1) we follow the recent approach of using fractional 



derivatives 



Id 



191 ] 20|, to generalise the Sturm-Liouville theory derivation of the GME in 



(2) and (3). This approach is based on the generalisation of the Brownian motion of the 
paths that are required for Euclidean Feymann path integrals in quantum mechanics to Levy 
random processes [2^. Rather than the GME being meromorphically continued into poles, 
in this case, the GME is meromorphically continued onto a annulus, corresponding to the 
entangled nuclear spin state. 0. Fonnany, tKe L.v,. patK .te,ra. ,e.e.tes a funct,„„a. 
measure in the space of either left or right continued functions, and has only discontinuities 
of the flrst kind (flrst order derivatives) corresponding to whichever spin-flip state of the ^^C 
nuclei is the most energetically favourable. We therefore replace the time derivatives in (2) 
and (3) by the following generalised (Caputo) fractional time derivative 2l|] 22|, 



Dl,, (5-(t)5+(0)) 



rri 



1 



d 



ly] Jo {t-t'Ydt 



- (5-(t')5+(0+)) 



(5) 



where S (t) and 'S'"''(0) are time evolution operators for the electron spin (in the adjoint 
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The relevance of this particular product 



representation) at time t = and t = t 
of evolution operators is that, following |6|] 9j 10|, the density matrix for the transverse 
electron spin in (1) is deflned via, Tr [ S'~(t)>S'^(0"'")S'^^(t)S'"^(0"^)] . Practically evaluating 
the Greens function propagators of this fractional derivative formalism then leads to overall 
constant in front the propagtors, compared to the conventional GME, since the Fourier 



transform of a fractional derivative is given by 24 1 



^ [D^, (5-(t)5+(0))] = -^^(t)^ [5-(t)5+(0)] , ^tit) 



V i sgn{t)8iT /2 



We 



(6) 



The fractional derivative and adjoint prescription we make formally deflnes a regulariza- 
tion of electron spin decoherence via a (^-function prescription 25 1. Our analysis of 



the decoherence of the electron spin polarization of Nitrogen Vacancy centers of Diamond ex- 
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tends previous analyses in two ways; firstly, we do not drop any off-diagonal (dipole-dipole) 



interaction terms in order to make the analysis more suitable for the GME approach in [13| , 
and secondly, we use the adjoint representation to explicitly calculate the nucleation of the 
of ^^C flip-flop dynamics via the factor l^l''. 



I. DISCRETE FOURIER TRANSFORM 



In the nonperturbative regime, the electron 



spin decoherence can be calculated for the 



GME in (1) using a range of numerical schemes 26|]; exact diagonalization [12] [27|[28|] , loop- 
cluster expansion |icl | and via a discrete Fourier transform. In the following analysis we 
use latter method (although the (^-function prescription can also be applied numerically via 



17( 1) since we are able to keep the long-time cutoff flxed in the discrete Fourier Transform 
via the lattice spacing, which allows us to make an accurate comparison of \t\'^ for different 
field strengths and polarizations. The advantage of this approach is in having in a proper 
regularization of the GME in a treatment which also includes a full dipole-dipole interaction 
and the field splitting of the Nitrogen Vacancy center. 

We calculate the decoherence of the transverse spin electron spin by using a discrete 
Fourier transform to calculate time-dependent Greens function of the Hamiltonian in (1), 
for up to a maximum of = 8 nuclei. Although = 8 is an experimentally relevant number 
for current entanglement studies 

000, this small bath size can further justified for larger 



26]. Exphcitly: if the 



quantum dot systems via the convergence of Chebyshev polynomials 
state vector of the spin-bath is represented via, 

Nb Nb 

\X) = J^a,\(3i/32... Pnb-1 ), 5^ |«fcr = 1 (7) 

k=l k=l 

where (3 is either or 1 (depending on whether the zth bath nuclei is spin up or spin 
down), ttfc is a random variable, and there are Nb = 2^ bath states, then, since the off- 
diagonal elements are of order 2~'^'^ / y/Ns, agreement with the electron spin polarization to 
6% of the value for the full system (of ~ 10^) can be obtained from sampling clusters of 
iust = 8 nuclei. Furthermore, this limit is an exact bound for the random Levy random 

n 

processes that form the (canonical) basis of our approach [29|. 

We obtain the energy-dependent Greens function of the Central spin model for Diamond 
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via an exact inversion of the Hamiltonian in (1) using Gauss- Jordan elimination, 



[G{e)],, = [{e - H)-%, (8) 

where G{e) is then a 2^^"''^ matrix, defined in the product basis = \S ^Ii^I2.--®In)- 
We repeat this calculation for M discrete values of energy, e, where (?„(£:) = G{e)\s=nAe with 
Ae = E/M and E the finite energy domain of the system, which we choose to be of ^(e^). 
The time-dependent Greens function of the Hamiltonian in (1) is then given by the following 
discrete Fourier transform, Gkit) = G(t)\t=k/\-r, which we evaluated using the numerically 



efficient Four-Step Fast Fourier Transform 30 1 



< r, g < p, < i?, g < P, pP = M. (10) 
where At = T/M with T the finite time domain of the system, and where the short 



range cutoff set via ET = h and M = 2 
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111 ]. The time-dependence of the transverse spin 



of the electon at the Nitrogen Vacancy centers in Diamond is then given by the expectation, 

{S+)t = {4GKt) S^G{t)ij,) = mm S-(t)5+(0+) IV^o) (11) 
where the initial electron state is defined to be -^n = ( |0) + |1) ) /a/2 perpendicular to the 

rinn 

initial direction of the nuclear spins, ^^M. Note that our choice of adjoint representation 
and fractional derivative approach means that the above expectation for the transverse 
electron spin is real (and negative) and is also not normalised to unity. This allows us to 
quantify the nucleation of the of ^^C fiip-fiop dynamics via the factor \tY in (6). 

We compare our numerical results in the first instance with the exact results for the 



Central Spin model obtained in 13| for the limits of large applied external magnetic field, 
and a fully polarized initial system. The zeroth order dynamics result for the transverse 
electron spin in (3) is given by, 

= (S+)oexp[-tV2r,2 + i{B{es - e,) +piV)t], r, = J (12) 

V^(i -r) 

where p is the fraction of initially polarized states, and Tc the transverse electron deco- 
herence time, which, from the perturbative expansion of (3), we can expect to be relevant 
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FIG. 1: Time dependence of the transverse component of electron spin, {S^)t, for spin baths of 
= 8 and = 6 ^^C atoms in Diamond, respectively. For N = 8 the applied field is varied 
between — lOOG, whereas for iV = 6 the applied field direction is varied relative to the z-axis at 
lOOG. 
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FIG. 2: Applied field strength dependence of the logarithm of the argument of the Caputo integral 
in (11), for spin baths of A^ = 2 — 8 ^^C atoms in Diamond. 

to order 2{k + 1) when the self-energy is supressed by a factor of {ujn/NY . From the second 
term in the exponent in ( fT2l) . it follows that the decoherence time, Tc, becomes independent 
of the applied field strength for large values of field, and that the nonperturbative regime is 
therefore delineated by, N\B-,\ > |^| ~ C(10 Gauss). 

In Figure 1 we plot our evaluation of the time dependence of the transverse spin com- 
ponent of the electron, {S^)t, for spin baths of = 8 and = 6 ^^C nuclei in Diamond. 
In the left plot (for = 8 nuclei) the applied field is varied between — 100 Gauss, and 
in the right plot (for = 6 nuclei at 100 Gauss) the angle between the applied field and 
2;-component of electron spin is varied between 9 = — n/2. Following our broad expec- 
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tations from the zeroth order result in (12), the decoherence time increases in the left plot 
from Ifis to 7/is as a function of decreasing applied field strength, and increases in the 
right plot as a function of increasing (polarization), 6, from 2/is to Qfis. Note, that whilst 
for N = 8 the expectation value of the transverse spin component of the electron tends to 
zero, (S'+)t^oo = 0, where we are closer to the above nonperturbative bound for = 6 the 
expectation tends to a nonzero value, {S+)t^oo 7^ 0. 

In Figure 2 we plot the applied field dependence of the logarithm of expectation value of 
the transverse spin component of the electron log((S'+)T) which from (6) gives a measure 
the relative differences in l^l'' as a function of applied field strength. From (4), near the 
upper bound of the memory kernel (z/ > 1) we expect a power law decay, where the electron 
spin decoherence envelope is slowly decaying and governed by the nuclei furthest from the 
central electron. As in Figure 1, evidently the systems for = 2 — 6 which are closer to the 
nonperturbative bound have larger associated values of u, which is also relatively insensitive 
to applied field strength. This is in contrast with the A^ = 8 system that is furthest from 
the nonperturbative bound. This implies that the meta-stability of the ^^C nuclei (which is 
responsible for the value of u) decreases as a function of increasing applied field for larger 
bath sizes, as would be expected from (12). 



II. SUMMARY 



We have calculated the transverse electron spin decoherence for Nitrogen Vacancy 
center in Diamond, modelled as a Central-Spin electron spin in the adjoint representation 
coupled to a bath of up to eight entangled ^^C nuclei, with an additional dipole-dipole 
interaction also included between these bath nuclei to enable us to quantify the flip-flop 
dynamics of the bath. Our approach is based on a (^-function regularization of the quantum 
tunnelling bound between the triplet and singlet states of the Nitrogen Vacancy center. 



2l| 



Our approach allows 



which we have introduced via fractional time derivatives 
us to extend previous analyses in the nonperturbative regime, which we have probed via 
the numerical Four-Step Fast Fourier Transform by calculating the decoherence of the 
transverse electron spin in (1) via a discrete evaluation of the time-dependent Greens 
function of the Hamiltonian. Our results for the decoherence time of the transverse electron 
spin are consistent with recent experimental measurements for (proximal) spin baths of 
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this size [ll-Qlj?], and furthermore, our investigation of the infrared cutoff dependence 
and regularization of our approach indicates that the stabihty of our analysis improves 
as a function of increasing bath size, as more channels become available for the quantum 
tunnelling. 

Our thanks to V.I. Fal'ko for useful discussions. 
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